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ABSTRACT 


Recent findings indicate that magnetic anomalies occur in rocks 
in the sector in which an earthquake is recorded. The enhancement 
of the local geomagnetic field is believed to result from the pressure 
buildup beneath the ground before the quake. Also, there are certain 
cases in which magnetic anomalies appear to correlate with areas 
of tectonic activity and mineral concentration. 

The antiferromagnet Haematite (a-Fe 2 0 3 ) exhibits an anomalistic 
behavior about a magnetic ordering temperature called the Morin phase 
transition. Below this temperature the system exists entirely in the 
antiferromagnetic state, while above it, a ferromagnetic phase, and 
hence a net magnetic moment appears. In this report, a more exten- 
sive spin wave theory of Haematite is developed. The behaviour of 
thermodynamic quantities around the Morin transition temperature is 
studied, and also the latent heat of the Morin transition is calculated. 
The temperature dependence of the antiferromagnetic resonance fre- 
quency and the critical fields Hc|| and Hci is calculated. It is 
found that the theory agrees well with experiment. 
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ON THE MAGNON INTERACTION IN HAEMATITE. 
II. MAGNON ENERGY OE THE ACOUSTICAL MODE 
AND MAGNETIC CRITICAL FIELDS 


I. INTRODUCTION 

All processes in nature including geophysical phenomena are thermodynamic 
in origin. Because of the magnitude of the pressures, volumes and tempera- 
tures associated with processes in the Earth's interior, most external mani- 
festations of such events are largo scale in time and space. 

Magnetic field studies are a source of information on the properties, structure, 
and dynamics of the Earth's crust, mantle and core. Geomagnetic anomalies 
are known to correlate with geologic structures and with areas of mineral 
concentrations and tectonic activity. One of the most important minei’al oxides 
to be found in sedimentary rocks is the magnetic compound known as Haematite, 
Haematite has the corundum structure and orders magnetically at the Neel 
temperature T N - 947 K. For temperatures between T N and the Morin tempera- 
ture Tm = 2G1 I<, haematite is a canted antiferromagnet with the spins perpen- 
dicular to the c axis. For temperatures below T^ , haematite is a uniaxial 
antiferromagnet with the spins lying along the c axis. The transition at T M is 
called the Morin transition. There are a number of experimental and theoretical 
studies most of which are described in the paper of Jacobs et al. (1971). 

The mechanism of the Morin transition in haematite was believed to be the sign 
change of the uniaxial anisotropy energy at T M (Kanamori 1963, Artman et al 
1965), However, several experiments including the antiferromagnetic resonance 
experiment of Foner and Williamson (1965) and the magnetically induced and 
the stress-induced Morin transition experiments (Flanders 1969, Allen 1973) 
were not explained by the simple anisotropy theory. In order to investigate 
the mechanism of the Morin transition, several people have recently studied 
the effect of the magnon interaction on the optical-mode magnon energy (Herbert 
1970, Nagai et al. 1972). On the other hand, the temperature dependence of 
the acoustical-mode energy was studied by Bonavito (1970) and Watarai (1973) 
by using the two-sublattice ractfsl for this crystal. Very recently, the present- 
authors (Nagai et al. 1973) 5 ‘.rm>isevv studied the spin-wave theory of haematite 
and pointed out that the critics 1 .perature T L at which the magnon energy in 
low-temperature phase becomes imaginary is slightly higher than the tempera- 
ture Tji at which the magnon energy in the high-temperature phase becomes 
imaginary. Thus, we have two kinds of locally stable states in a narrow tem- 
perature region near T M . 
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II. ACOUSTIC MODE ENERGY 


We start from the following model Hamiltonian 


11 = 2 J £ S,.S„, 
o .«•> 



D Y. ( SjxS m U 


assuming the two-sublattiee body-centred lattice. Here Sj and S m are the 
vector spin operators associated with the jth and mth atoms of the crystal and 
the sum i^ m .is extended over all nearest-neighbour pairs. The exchange 
interaction is given by the first term of (1). Since we assume J > 0, the jth 
and mth spins are coupled antiferromagnetically. The second term denotes 
the one-ion type uniaxial anisotropy energy which makes spins direct to tl < * / 
axis. On the other hand, the dipolar type interaction which is given by the 
third term tends to make the spins lie in the xy plane. The role of the 
Dzyaloshinsky-Moriya (DM) interaction, the last term, is similar to that of 
the dipolar-type interaction except that the DM Interaction makes the interact- 
ing spins orthogonal to each other when the spins lie in the xy plane. 

Recently, Jacobs et al. (1972) pointed out that the crystalline anisotropy energy 
of the form 


- C[£ S% + £ 2 

may be important in explaining AFMR experiments. The effect of this energy 
will be discussed in rv. 

In the low-temperature phase, the spins are parallel to the iz axis, as described 
in the Introduction. By using the method described in our previous paper (1973), 
the magnon Hamiltonian II L is given by 


H h = (2JSz)Y hMok + W + gzViflkbk + oivifllbl + gtfMb-k + 

k (3; 

+ 0A a k a -k + bk^-k) 4 * + blbLk)], 
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Hero g, etc are given by 

= (1 - — 'V Lf ) - e(l — »l) + v 0 " 1/25 — 2iii) - 'uiw u , 

0 z = (1 - »l ~ wj - if/d - Hi) + Ov L% (4) 

y 3 = ~x L (l - e) - JO'l + y£) - M/l - >’l). 

04 - ~{{2vy L + Xi - ulx l). 

where e = A/2J and i> - B/Jz. The parameters u L , w L , x L and y L are defined 
by (2/NS)S <at,a q >, (2/NS)S q < a^bf, > 7«,, (2/NS)S q < a + q b H ,>y q and (2/NS) 

S q <a q a. q > respectively. In (4), w u and w Li are defined by (w L + w+)/2 and 
( Wl - %v^)/2, respectively. The diagonalization of the Hamiltonian may be done 
as follows. Confining our attention to the positive half space of k, we define 
tho new operators: 


« +1 = (S 1,2 (^» + c 0. b-jt = - O. 

- & m {B k + DO, 6 t , = - DO. 


( 5 ) 


Furthermore, we introduce F k and G k by 


, cA k + 2gt (J 2 k - 2 c M 

A k -\- — r^n J b w - “ 


( j2k + 2(74 | 


1 02k “* 2 04 


where g 2k = g 2<rk and g 3k = g 37k . 

Finally, introducing rhe operators a lk and 0 lk by 

Ft = aU cosh 0 k - /i lk sinh 0 k , B\ = /?U cosh 0 k - a,* sinh 0 fc , 
GJ = a J k cosh - (i 2k sinh </>*, D[ = /i£ k cosh - a 2k sinh </> k , 


where 


tanh 20, 


= and tanh 20* = ^ , 


+ 03k 


01 “ 03k 
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wo obtain 


J ( Z tmMk*n + i), 


( 0 ) 


where 


/k% = !/&[(*, - (-I )'</ 3A ) 2 - \g u - (-1 )' 20 4 |*] l « 


(7) 


In our approximation, u L etc are determined from the following equations: 


" L = 2NS ^ ^ C ° Sh 2 ° lk{l + 2/ '^ + COsh 2 ° 2k(l + 2n ^ - 2 l 


U ’ L 2 NS? 


9*k -1- 2gj 

. k 


x (1 + 2n\J + — M x (1 + 2 /jjJ 

' l(tJ 2k 


-_L_v 

Al 2 NS? 

1 r 
3 ’ l 2 NS ? 


f/i + 173* 


tuo lk 

02 ft + 2$f* 


X (* + 2 |«V*) - X (1+ 2/i 2 J 

fta) 2k 


x (2 JSz), 
( 8 ) 

x (2,/Sz), 


ha) 


i* 


x (i + 2**0 - 2JL--M X (1 + 2.1*0 
hcv 2k 


x (2 JSz), 


where n^ - [exp (/3hoj ik ) - 1 ]' J with 0 - l/k u T. In the magnon renormaliza- 
tion approximation (MBA), u L , w L , x L and y L are calculated from (8) self- 
consistently. 

Since u L etc are of the order (1/S), the square root of (7) can be expanded in 
a power series with respect to (1/S). Since we retained only those terms up 
to order of (1/S)* 1 and (1/S)° in the Hamiltonian, the square root of (7) will be 
accurate up to order of (1/S), as pointed out previously (Nagai 1969). 


By assuming that the anisotropy energy is very much smaller than the exchange 
energy, the temperature-dependent magnon energy with wavevector k can be 
written as 


4 


hu) ik = 2JSz[(\ - yl){\ - u L - w L ,) 2 + (1 - - w lf ) {2v'(l - 2»J (9) 

- (2e + d 2 ){ 1 - ti L + Wl,)}] 1 ' 2 , 

whore v 1 ~ v(\ - 1/2S). For the sake of convenience, we give the estimation of 
the magnitude of the anisotropy parameters in (9). According to the experiment, 
the critical field II c |j *s 68 IcG at T -- 0 K. This corresponds to htojo - 10 K. 

We may assume 2JSz 1000 K because T N ■- 947 K. Hence wo find that 
(2 v 1 - 2e - d 2 ) is of the order 10" 1 . On the other hand, the experimental value 
of W c i is 160 kG at T - 0 K, The magnitude of d is estimated to be of the order 
of 10' 2 from the relation d ~ (gp B /2JSz) (II 2 ||/iI c i), which will be explained in 
III-2. The calculation of dipolar energy leads to e =■ 10’ 3 , and hence v is sup- 
posed to be of the order 10 * 3 . 

If we use (9) for the magnon enei’gy, the equation for T L at which ho> i0 (T) vanishes 
is given by 


2v'(l - 2 m J - (2e + f / 2 ) ( 1 - m l + wj = 0. 


By using this equation, T L is calculated as 


(k a r L \ _ (l5S\ l, *(2v' - 2e - d 2 \ u * 

\2fSzJ~\4n 2 ') \ 2e + d 2 ) ’ (U) 


and hco 10 (T) is written as 


luo i0 (T) = hw LQ (T) = hco L0 { 0) [1 - (r/TJ 4 ] (12) 


where hco LO (0) = 2JSz(2y' - 2e - d 2 )‘ /j . 

In the high-temperature phase, the spins are in the xy plane and they are 
canted slightly by the DM interaction. If we denote the angle between the jth 
and the mth spin vectors by (n -2 0) as seen in Figure 1, the Hamiltonian is 
written as 


5 


X 


Figure 1. 

//= I [(2 J cos 20 + D sin 20) + S’„6' nl „) + (27 - ,1) 

O.m) 

+ (27 sin 20 - D cos 27) [S A S m - S jn S )] - Dfc S} 4 + £ S&. 

J m 

(13) 

The angle Q may be determined from the condition 2J sin 20 - D cos 20-0. By 
using again the Holstein-Primakoff formalism, the magnon Hamiltonian for the 
high-temperature spin configuration is given in the same form with (3). In this 
case, g 3 etc, which are written as g m etc, are given by 



i 

i 


On t — (1 — »n — w Jlf — iv,j) 4- ic(w Hf + ,v H ) + \(l~{ 1 — U|i — H’ IIr + j,v fl ) 
- £v*[l - 2 «„ - &y u + rfi)], 

On 2 - 0 - »n - "'ll) — i«U - tin “ W + W 2 [iU “ »u) - w„ + J >•„], 

0n5 - -- x n - iOu + y,v - M 1 - »h ~ u(.v» + jriM - i^ 2 [i(i - %) 

+ Xu + aO’lt •+■ >■*[)]. 


On* — — 8 (2.V,, + U’li) + + W,j) + i t (l 2 (W u — .V lt ) — 


h‘V + 1/4S - in,, - i y lt ) 


t 

f 


I 
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Wo may calculate the magnon onor&y hou^T) using tho formula (7), and then 
wo obtain the critical temperature T|| whore tho magnon energy with infinito 
wavelength vanishes, hw, I0 (T Jt ) - 0. Wo havo shown in tho previous paper 
(Nagai ot al. 1973) that the critical temperature in the low-temperature phase 
T l is slightly higher than tho critical temperature hi the high- temperature 
phase Tjj. Tho difference botweon them is obtained as x T L and it is of 
the order 0* 1 K. Using tho condition that tho free energies in both phases aro 
equal, wo calculate the Morin temperature T M as T M - T L x (1 - £n). By using 
tills T m in (12), hejo (Tm) is estimated to be only a few kG, According to Fonor 
and Williamson (19G5), tho magnitude of hwo(T M ) is 28 kG. Thus the theory 
does not agree with tho experiment on this point. Recently, Jacobs et at. (1971) 
pointed out that this contradiction may be removed if we assume the S, 4 type 
crystalline anisotropy energy. Evon if this type of anisotropy is introduced the 
magnon energy spectrum with k 4 0 may not be affected and accordingly the 
thermodynamic quantities are not affected. Hence we calculate the thermodynamic 
quantities without introducing the higher-order anisotropy enorgy in this section. 
The effect of this anisotropy energy on the physical quantities will be discussed 
in IV. 

It was shown in our previous paper (1973) that the magnon energy with k 0 at 
moderately high temperatures is given by 


haUT) = H u (l - y/)'» (15) 

where H L 2 - 2JSz(l - Ul~ Wl.) for the low- temperature spin configuration and 
H U 2 ~ 2JSz(l — u j j — w,j - £e) for the high-temperature spin configuration. The 
thermodynamic quantities around the Morin temperature can be calculated in 
the same way as in a previous paper (Nagai and Tanaka 1969). The difference 
between the sublattice magnetization M just below T M and M just above T M is 
calculated as 

/ /• r P \ 2 

AM = M(T m - 0) - «rr a + 0) a -J X v. m 


The corresponding quantity in the case of the magnetic specific heat C(T) is 
calculated as 

j 

ti 

n 


7 


/o n 2\ f U r r \3 

ac * qr„ - o) - c(7- m + o) = Nk a 


X V. 


(17) 


Thus, wo find both (AM/M) and ( AC/C) arc smaller than 1 ft. Recently, Allen 
(1973) analysed the experiment of the stress-induced Morin transition of this 
substance. He estimated the latent heat of this transition to be 5 x 10 s erg/cm' 3 . 
According to tho present theory, the latent heat is calculated as 


7;, x [S(7 m + in ~ S(T sl - 0)] - A%7 M 



(18) 


where S(T) denotes tho entropy. The latent heat of (18) is evaluated to bo 
2 x 10 8 x v erg cm' 3 . Tills value is of the same order of magnitude as tho 
Allen value. In the above calculations, we neglected the effect of the optical- 
mode magnon energy because its magnitude Is of tho order 10 3 K. 


III. MAGNETIC CRITICAL FIELDS 
1. Parallel Critical Field II c j| (T) 

It is known that tho spin-flop critical field H c j| in the case of the uniaxial 
antiferromagnet is, in the molecular-field approximation (Nagamiya et al. 
1955), given by 


Wtnllc n 


= hio 0 (T) 



(19) 


On tho other hand, the spin -wave theory gives the same result at low 
temperatures (Nngai and Tanaka 1970), In the latter paper, the critical 
field is defined as the field at which the energy of the magnon with Infinite 
wavelength vanishes. Exactly speaking, the critical field II C || obtained in 
this way is slightly diffei’ent from the critical field obtained from the con- 
dition that the free energies in both the flopped and unflopped phases are 


S 


oqual. Actually, according to our theory, APMR does not vanish at T M and 
Its magnitude is calculated to be 2 kG at T M . Therefore, Ii c || (T M ) is also 
equal to 2 kG if the equation (19) is used for II C || . However, it is known 
that II C || (T) must vanish at tho Morin transition temperature because the 
Morin transition occurs without a magnetic field. If we inspect the equa- 
tion for determining II c j| (T) more precisely, wo will find II C || (Tm) = 0. 
Actually, T L is very close to T M and !I C || (T M ) = 2 kG is very small com- 
pared with II C || (0) = G8 kG, Hence II C jj (T) is approximately given by (19) 
of which tho main temperature behaviour is written as 


M e »(T) = H cll (0)[l - (T/T,,) 4 ] 1 ' 1 (20) 


whore II 0 || (0) = (2JSz/gju B ) (2v ’ - 2e - d 2 ) w . The theoretical result of 
II C || (T) is compared with experiment of Foner and Shapira (1969) in Figure 2. 



Figure 2. Temperature Dependence of H c || (T) which is Computed from 
(19) is Compared with Experiment (Foner and Shapira 1969). 
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Next let us calculate and H tf || (T) in the njulecular-fiold approximation 
In order to seo the correspondence of the spln-wavo result to the molecular- 
field rosult. Wo neglect the DM interaction for simplicity, 

The Weiss molecular-field Hamiltonian IljJ in the unflopped phase Is given 
by 


«m ^(2J-A) Y [S J( <S„ ( >+ <S Jt > S„. - < S, ( ><S„ ( >] 

<;.m> 

- B [£ S ,1 - r £ si.] - mil [£ S, ( + £ S„ ( ], (21) 


Let us introduce the variables p and n by < S Jt -> - p + n and <S II){ .> -■ 

-p + n. Here p denotes the average value of jth spin in the cfcso If = 0 and 
n is approximately proportional to If, Furthermore, wo define p 0 and p' 
by P = Po + p'. Here p' is proportional to A or B and p 0 - <Sjj.> 0 -• 

Tr tSjj-exp{- /JHm ) ] /Zo» where H®j = -2Jzp 0 Sj Sj f and Zo ~ Tr [ exp(- (3 Hj(j ) . 
free energy Fjtj in the unflopped phase Is calculated as 


= -k u T\n Z 0 + NJzpK 1 + e) - A7rv<5j> 0 - (% ( ,// l /2). (22) 

The molecular-field Hamiltonian II m in the flopped phase is given by 
H ( m = £ C(2 J cos 20 + A sin 2 0) {S, ; <S > + <S > S - <S, ( ><S >} 

- (£/r){(S 2 . + Sty sin 2 0 + (S} { + S 2 i4 ) cos 2 0 }] (23) 

j m 

where '.he angle between II and Sj is given by (&7T-0). 

Let us define p and n' by < Sjj.> = -<S mf > = p + n' where p is a 0 -indepen- 
dent tern. Again we introduce p 0 and p" by p = p 0 + p" where p" is pro- 
portional to A or B. Differentiating the free energy f/j with respect to 0 , 
we obtain sin 0 = gpu Il/4Jzp 0 . Finally, the free energy is written as 


f f M = NJzpl - k ti T In Z 0 - - (y. x H 2 /2). 


(24) 



By equating F^ to , we may calculate the Morin temperature T M and 
also the critical field II C || (T). 

The equation for determining T M in the case II = 0 is given by 


2 v<Sj f - £%> o - 2 €<S J; ) 2 0 = 0. 


(25) 


This equation was already used previously in order to explain the Morin 
transition (Kanamori I9f?3). The t itical field II C || (T) is calculated from 
the following equation: 

«- 0 H tll (T) = 2J:rjv<S* x - Sj<>„ - &<s*>sr x [I + fe,/2xJ]. ( 20 ) 

wliicn should be compared with (19). If the DM interaction is taken into 
account, e in (25) and (26) is replaced by fe + (d 2 /2) ] . 

2. Perpendicular Critical Field H c i 

If a magnetic field is applied perpendicular to the c axis, say, parallel 
to the y axis, then the spins rotate in the xz plane and suddenly flop into 
the x axis at the critical field H c x (T). This situation can be seen in 
Figure 3. 



Figure 3. Spin Arrangement in the Magnetic Field (H|ly Axis). 
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The Hamiltonian of the system is tne sum of (1) and the Zeeman term 


l 

| 

j 

i 

! 

f 


Wz “ - rn » H (Z s Jy + Z s m>)‘ (27) 

i m 


We take the coordinate axes |£jT?j f jj for the jth spin and the axes j £ m j 
for the mth spin. In this new system Sj and S 1U are written as 




4MM) 



(28) 


whore 


0±(fl, 0) 


( ,COS (f) 

.0 

-sin (p 


+ sin Osin <p 
cos 0 

+ sin 0 cos (p 


ms 0 sin <fi 
±sin 0 
COS 0 COS (j) 



(29) 


Substituting (28) into the Hamiltonian and expanding the Hamiltonian in 
terms of spin-wave operators, then we obtain a series 


H - E 0 + H x + H i + II 2 + II 4, (30) 


where E 0 denotes the classical energy and is of the order S 2 , and H lf II 2 , 
II 3 and II4 are, respectively, of the order S 3 ^ 2 , S 2 ^ 2 , S I/,2 and S°. 

If we omit the last two terms in (30), we can obtain the tempei-ature- 
independent solution of H c _l . The angles 0 and 0 are, in this approxi- 
mation, given as follows; (i) sin 6 = h (u - e )/(2v - 2e - d 2 ) and sin 0 = 
d tan 6/( v- e) (unflopped phase); (ii) sin 0 - (h + d)/2 and 0 = n/2 (flopped 
phase). I-Iere h is defined by h = gju B H/2JSz. 

The free-magnon Hamiltonian in the unflopped phase is obtained in the name 
manner as that given in the previous section. After a somewhat extended 
manipulation, the free-magnon energies are given by 
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Uo)\ k = 2JSz\\ - y\ + 2v - 2e - d 2 - (d 2 h 2 /(2v - 2c - c/ 2 ))] 1 /2 , 

* = 2JSs[l — y 2 4- 2v — 2c — d 2 + /i 2 {4(v — c) 2 — 4f/ 2 (v — c) 

-</ 4 }/(2v ~2c -d 2 ) 2 ]' 12 , 

where we used the relation 1 > d 2 , v , e . The critical field h c i is deter- 
mined from the condition luoio = 0 and is given by h c i = ( 2 r - 2 e - d 2 )/d. 
On the other hand, h c j| is given by h c || = (2 u - 2e - d 2 ) 1/1 and hence h c j. is 
written as h c i = h c 2 || /d. 

The magnon Hamiltonian in the flopped phase is obtained similarly. The 
free-magnon frequencies in the flopped phase are written as 


hco[ k = 2JSz[[ — y 2 + hd — (2v — 2c — t/ 2 )] 1/2 , 
ha) ( 2k = 2JSz[l - y 2 + h(h 4* i/)] w2 - 


The temperature dependence of the critical field H c x (T) may be obtained 
from the condition that the temperature-dependent magnon energy with 
inifinite wavelength vanishes. The calculation of H c j. (T) was performed 
as follows: (i) the temperature dependence of 0 and 0 is taken into account 
by applying the RPA in the cubic terms of magnon operators in H 3 ; (ii) the 
quartic terms of magnon operators in H 4 are treated by using RPA; (iii) the 
non-diagonal matrix elements of the cubic terms of magnqn operators in 
H 3 are taken into account by using the second-order perturbation approxi- 
mation (Maekawa 1973). After a complicated manipulation, the tempera- 
ture dependence of the perpendicular critical field was found to be h c l (T) = 
[ h c || (T) ] 2 /d. We will not writ e out the derivation of this result by the use 
of the above-mentioned method, because it is much too complicated and 
lengthy. Instead, we derive the same result by the use of molecular-field 
approximation. 

The molecular-field Hamiltonian is written as 

H m -2J £ [(cos 29 - cos 2 9 cos 2 <j) + d sin 29 sin 0) {S } ^p m 4 - - PjP m ) 

- (e/2) {sin 2 0(SJ 4- S 2 J + sin 2 0 cos 2 <t>(S% + S 2 ^} (33) 

-ShrinOiSjt^SJl 

where Pj = < Sj r >. By using a method similar to that given in the previous 
section, the free energy in the molecular-field theory can be calculated. 
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From the conditions 3F M /30 = 0 and 3F M /30 = 0, we obtain the following 
two equations: 

cos 0 cos $[cos f? sin 0(v<S^ — S^) 0 - e<S, ; >n) — d sin 0<Sj f >o] = 0, 

sin 20[2<S j{ >o + cos 2 0{v<S^ - S 2 { ) 0 - €<S, ; >q}] -2d cos 20 sin <KS J{ > l (34) 

= 2S<Sj ; ) 0 h cos 0. 

In deriving (34), we used the relation <Sj f + Sj„ Sj t > 0 = <Sj nt S, m) + S„ ln S me >o 
= 0. In the unilopped phase, 0 and <fi are determined from the following equa- 
tions: 

sin (J) - d tan 0j{\\ - e), 

h = sin 0«S J( >o/S) x (l - (35) 

ty 

where u l - v<S# - Sj ? ^/< Sj f >o . On the other hand, we have </> = n/2 
and 


2<S J{ > 0 sin 0 = Sh + rf<S ;{ > 0 

(do) 

for the flopped phase. By equating F M (flopped phase) to Fm (unflopped phase), 
the equation for determing the perpendicular critical field is given by 


h' = (2 Vi - 2c - d 2 )/d (3 7 ) 

where h* - Sh/<Sj f > 0 . If we use (26) we may write 

*‘ 1<n = 5<s^ [A ‘ l|(:0]! < 38 ) 

of which the main temperature dependence is given by [1 - (T/Tm) 4 1 . The 
computed result of h c j (T) is compared with experiment in Figure 4. 


IV. FOUR- SUBLATTICE MODEL 

Let us compare the theory based on the four-sublattice model with experiment. 
As can be seen in the previous paper (Nagai et al. 1972), the energy of the 
magnon with infinite wavelength is written as 

(Iiioq) 2 = - S 2 D 2 + 2S 2 (2J 0 + 6J 2 + 12J 3 ) [2 V - ^0> u - <D I2 - <1> U + <h 14 )]. 

(39) 
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Figure 4. Temperature Dependence of I-I 0 j. (T) which is Computed from 
(38) is Compared with Experiment (Foner and Shapira 1969). 

where D denotes the DM constant and U the uniaxial anisotropy constant. The 
dipolar anisotropy constants 0 NM are defined by 

£ (RJ„ - 3 Z^,)/(R„,) S = (4/i S/a 3 ) P m 

«NM 

where a denotes the lattice parameter 5 • 42 x lO'^cm, Rn and Rm belong to 
the N and M sublattices, respectively, and Rnm = Rn “ e m- The numerical 
values are given by P^ - 1 - 09, P ]2 =*-9*15, P 13 = -2 • 03 and P l4 = 11 * 01. By 
using the exchange parameters proposed by Samuelsen and Shirane (1970), we 
obtain 2J 0 + 6J 2 + 12J 3 = 468 IC. It is seen by inspection that we may adopt 
the following substitutions: 

2 Jz — 2 Jq 4- 6J 2 4* h — (uPhH/S)/{2Jq 4~ 6J 2 4* 12 ^ 3 ), 

v = U/(J 0 4- 3 J 2 4- 6 J 3 ), d = D/{2J 0 4- 6J 2 4- 12 JJ, 

— 3 ( l ) n — O 12 ~~ ^13 4- O 14 
€ ~ 4 2J 0 + 6J 2 4- 12 J 3 
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By using the experimental values of the critical fields at T = 0 K, wo obtain 
U = 0 • 222 k u and T M is estimated to bo 258 K. 

1. Effect of Sj-term 

In tills section, we give a qualitative discussion of the effect of the higher- 
order anisotropy energy on the magnon frequency. If wo use the Holsteln- 
Primakoff formalism for spin operators, the higher-order anisotropy 
energy is written in a series 

-CS 4 = -CS 4 [ 1 - (4/S) (1 - 3/2S) a]cij + (6/S 2 ) a^aj 
-f higher-order terms of (1/S)]. 

We retain up to the thix’d tei'm in the brackets and neglect the others. By 
using the method described in II, the magnon energy of (9) can be written 
as 


hcoUT) - 2JSz{(l - y k 2 )(l - u L - w L ) 2 + (1 - u L - w L ){2v'(l - 2ud 

+ 4c(l - 6w l ) - (2e + d 2 )(l - u L -f vv L )}] 1/2 , ( 40 ) 


where c = CS 2 (1 - 3/hS)/Jz. The critical temperature T£ at which hw L0 (T) 
vanishes is given by 


k B T' L \ A 5S\ 1 ' 4 f 2v' -f 4c — 2e — d 2 \ 1 /4 

2JSz)~ W 2 / \ 2 € + d 2 ) 


On the other hand, it is easily seen that the magnon energy in the high- 
temperature phase is given by 

/ico, u (T) = 2JSz[(I — }'*)(1 — u H — w„ — ^) 2 + (1 — Uj| — vv tl ){(2€ + t/ 2 )(l — w„ + vv,,) 
— 2v*( 1 - 2 «h) - 3c(l/25) + «h)}] ,/2 . (42) 


Thus, the magnon energy hco Ilk (T) is not affected by the introduction of c, 
and it is approximately written as 

/ia> H0 = 2JSz(2v' ~2t- d 2 ) U2 [- 1 + (T/T„) 4 ] l/2 , 

(43) 
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where 


fk 0 T u \_ p5SY /4 /2^-2€-rf 2 Y M 

V V 2€ + rf 2 ) ’ 

The Morin temperature can be calculated in a manner similar to that used 
in ni-1. The equation for determining T M is written as 

2v'(l - 2ui) — (2c + d 2 ) (1 — w L + wj + 2c(l - 4 »l) <= 0, (44) 


by which T M is written as 

/ WA_/15SV'Y 2/ + 2c - 2e - </ 2 
\2 JSz) \4a 2 J \ 2 € + d 2 

If we ULe the above expression for T. ( , AFMR is written as r 

p 

|i 

ho) LQ (T) = 2JSz[2c + (2v* + 2c - 2e - d 2 ){ 1 - (T/T N1 ) 4 }] W2 , (46) Ij 

li 

11 

v f [ 

and we also obtain ho>Lo(T M ) = hw H0 (T M ) = 2JSz(2c) ”, By using the experi- jj 

mental value hto 0 (T M ) = 28 lcG, c is estimated to be 0 ‘ 58 x 1Q' S . If we f[ 

write (46) in the following form 


hco L0 (T) - /im LO (0)[I - (T/TW 2 , (46) 

T 'l is written as T^ = T M (1 + 5 ). Here 5 = (c/2)/(2 v' + 4c - 2 c - d 2 ) is 
estimated to be 0 • 04. The temperature dependence of hca L0 (T) which is 
computed from (40) is compared with experiment in Figure 5. Thus, the 
critical temperature T^ and also hco H0 (T M ) are strongly modified by the 
introduction of c. However, the magnon energy spectrum does not strongly 
depend on the magnitude of c at moderately high temperatures, as can be 
seen from (40) and (42), Therefore, the temperature behaviour of thermo- 
dynamic quantities may not be affected by c. The previous result (16) ~ (18) 
will still hold. Here T M in these equations should be replaced by T M of (45). 
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tck) r M 

Figure 5. Temperature dependence of hcojjo(T) and hoj| JO (T), which are com- 
puted from (40) and (42), respectively. Open circles are deduced 
from 70 GHz data (Foner and Williamson 1965). 


Similarly, the critical fields li c || (T) and h c x (T) are, respectively, given 
by (20) and (38), where h c |( (0) and T M should be replaced by (2o’ + 2c - 
2e - d 2 )’ /j and (45), respectively. 


\ 

i. 

* 


V. CONCLUSION 

We have seen in this paper that the calculated result of H c jj , H c x and hw 0 (T) 
shows a very good agreement with the experimental result. We took not 
2JSzt2u (1 - 1/2 S - 2u l ) - (2e + d 2 ) (1 - u L + w L )] l/j but 2JSz[2t>(l - 1/2S) (1 - 2u L ) 
- (2e + d 2 ) (1 -u L + w L )] Vi for the expression of hco L0 (T). These two expressions 
are both accurate because the square root is accurate up to order of (1/S), as 
noted in n. The latter expression shows [1 - (T/T m ) 4 ] Vi dependence, while 
the former shows a slightly different temperature dependence. The S 2 term 
does not work as the anisotropy energy in the case S - and hence the latter 
expression will be more adequate. The derivation of the latter expression for 
hco L0 (T) from the theoretical point of view may be made by using a more 
accurate theory than the present RPA. This remains to be done as a future 
problem. 
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